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CONVERSION  TABLE 

Conversion  Factors  for  U.S.  Customary  to  metric  (SI)  units  of  measurement. 


MULTIPLY  - ^  BY  - ►  TO  GET 

TO  GET  ◄ -  BY  ◄ -  DIVIDE 


angstrom 

1.000  000  X  E  -10 

meters  (m) 

atmosphere  (normal) 

1.013  25  X  E  +2 

kilo  pascal  (kPa) 

bar 

1.000  000  X  E  +2 

kilo  pascal  (kPa) 

barn 

1.000  000  X  E  -28 

meter2  (m2) 

British  thermal  unit  (thermochemical) 

1.054  350  X  E  +3 

joule  (J) 

calorie  (thermochemical) 

4.184  000 

joule  (J) 

cal  (thermochemical) /cm2 

4.184  000  X  E  -2 

mega  joule/m2  (MJ/m2) 

curie 

3.700  000  X  E  -1 

*giga  becquerel  (GBq) 

degree  (angle) 

1.745  000  X  E  -2 

radian  (rad) 

degree  Fahrenheit 

tr«(t°f  +  459 . 67 ) /I . 8 

degree  kelvin  (K) 

electron  volt 

1.602  19  X  E  -19 

joule  (J) 

erg 

1.000  000  X  E  -7 

joule  (J) 

erg/second 

1.000  000  X  E  -7 

watt  (W) 

foot 

3.048  000  X  E  -  I 

meter  (m) 

foot -pound- force 

1.355  818 

joule  (J) 

gallon  (  U  S  liquid) 

3.785  412  X  E  -3 

meter3  (m3) 

inch 

2.540  000  X  E  -2 

meter  (m) 

jerk 

1.000  000  X  E  +9 

joule  (J) 

joule/kilogram  (J/kg)  radiation  dose 

absorbed 

1.000  000 

Gray  (GY) 

kilotons 

4.183 

terajoules 

kip  (1000  Ibf) 

4.448  222  X  E  +3 

newton  (N) 

kip/inch2  (ksi) 

6.894  757  X  E  +3 

kilo  pascal  (kPa) 

ktap 

1.000  000  X  E  +2 

newton-second/m2  (N-s/m2) 

micron 

1.000  000  X  E  -6 

meter  (m) 

mi  I 

2.540  000  X  E  -5 

meter  (m) 

mile ( international ) 

1.609  344  X  E  +3 

meter  (m) 

ounce 

4.448  222 

kilogram  (kg) 

pound-force  (lbs  avoirdupois) 

1.751  268  X  E  +2 

newton  (N) 

pound-force  inch 

6.894  757 

newton-meter  (N* **m) 

pound- force /inch 

2.834  952  X  E  -2 

newtodmeter  (Nlm) 

pound- f orce/foot2 

1.129  848  X  E  -  1 

kilo  pascal  (kPa) 

pound-force/inch2  (psi) 

4.788  026  X  E  -2 

kilo  pascal  (kPa) 

pound-mass  (Ibm  avoirdupois) 

4.535  924  X  E  -  1 

kilogram  (kg) 

pound-force/foot2  (moment  of  inertia) 

4.214  011  X  E  -2 

kilogram-meter2  (kg*m2) 

pound-mass /foot3 

1.601  846  X  E  +  I 

kilogram-meter3  (kg*m3) 

rad  (radiation  dose  absorbed) 

1.000  000  X  E  -2 

**Gray (Gy ) 

roentgen 

2.579  760  X  E  -4 

coulomb/kilogram  (C/kg) 

shake 

1.000  000  X  E  -8 

second  (s) 

slug 

1.459  390  X  E  +  I 

kilogram  (kg) 

torr  (mm  HG.  0°C) 

1.333  22  X  E  -1 

kilo  pascal  (kPa) 

*The  bacquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  =  1  event/s. 

**The  Gray  (GY)  is  the  SI  unit  of  absorbed  radiation. 
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SECTION  2 


DATA  COLLECTION 

Our  data  set  of  the  Annual  Bulletin  of  Chinese  Earthquakes  (ABCE)  consists  of  both  paper 
bulletins  and  digital  data  sets  provided  by  our  colleagues  from  the  China  Seismological  Bureau 
(CSB).  These  data  are  from  over  1 10  stations  of  the  Chinese  national  network  (Figure  2-1 ).  The 
seismometers  are  mostly  short-period  (~1  Hz)  recording  on  paper.  Long-period  instruments  with 
periods  of  around  1 0  s  are  usually  co-located  with  the  short-period  instruments.  These  were  used 
for  magnitude  measurements. 

We  obtained  paper  copies  of  the  ABCE  spanning  1983  to  1995  and  digital  versions  of  the  1991 
to  1995  bulletins.  The  digital  data  contain  the  same  locations  and  information  as  the  paper 
bulletins  but  often  include  readings  at  international  stations  that  were  collected  by  the 
International  Seismological  Centre.  We  had  data  from  the  1984  to  1989  bulletins  transcribed 
into  digital  form.  This  yields  over  ten  years  of  digital  data. 

We  also  collected  paper  copies  of  the  Bulletin  of  Seismological  Observations  of  Chinese  Stations 
(BSOCS)  for  1979,  1992,  and  1993.  These  bulletins  are  a  teleseismic  version  of  the  ABCE. 
They  contain  mostly  readings  from  about  25  Chinese  stations  that  traditionally  report  to  the 
International  Seismological  Centre.  Most  of  the  readings  are  teleseismic,  but  local  events  are 
also  included.  Event  locations  are  either  from  the  ABCE  or  the  ISC. 

The  data  quality  from  the  ABCE  is  exceptional.  Seismic  stations  in  China  have  routinely 
identified  and  picked  all  seismic  phases  both  regional  and  teleseismic  for  events  in  and  around 
China  and  used  them  for  event  location.  These  phases  include  Pg,  Pn,  P,  PP,  pP,  Sg,  Sn,  S,  SS, 
sS,  sP,  PcP,  SI  1,  ScS,  PI  1,  ScP,  and  PcS  (Figure  2-2).  Regional  Pg,  Pn,  Sg,  and  Sn  phases  were 
identified  up  to  nine  degrees  distance,  after  which  they  were  classified  as  S  or  P.  In  addition  to 
travel-time  picks  and  locations,  the  amplitude  and  period  measurements  needed  for  mb,  mB,  and 
Ms  magnitude  estimates  were  also  routinely  reported  by  the  ABCE.  Amplitude  measurements 
were  routinely  made  at  the  majority  of  stations.  The  quantity  and  thoroughness  of  these  bulletin 
data  make  this  the  best  such  data  set  in  the  world. 
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Figure  2-1.  Shaded  topography  of  China  showing  locations  of  regional  seismic  stations  used  in  the  Annual  Bulletin 
of  Chinese  Earthquakes. 


Figure  2-2.  P-wave  first  arrival  travel  times  from  the  Chinese  bulletins  (ABCE)  for  1991  through  1995.  Pg,  Pn  and 
a  distant  P  branch  are  clearly  visible.  The  ABCE  routinely  identifies  arrivals  as  Pg,  ,  P,  PP,  pP,  Sg, 
Sn,  S,  SS,  sS,  sP,  PcP,  SI  1,  ScS,  Pit,  ScP,  and  PcS. 
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SECTION  3 


PN  TOMOGRAPHY  BENEATH  CHINA 


3.1  INTRODUCTION. 

We  have  applied  the  Pn  tomography  algorithm  of  Hearn  (1999)  to  Pn  data  from  the  ABCE  to 
image  the  P-wave  velocity  and  anisotropy  at  the  surface  of  the  mantle.  Raypaths  for  these  data 
cover  most  of  central  and  Eastern  China.  Tibet  and  Western  China  have  scant,  but  useable 
coverage.  We  used  raypaths  between  1.8  and  15  degrees  distance  and  only  stations  and  events 
with  more  than  9  arrivals  at  them.  The  final  Pn  data  set  used  contains  over  50,000  arrivals  that 
span  most  of  China  (Figure  3-1). 


Figure  3-1.  Pn  raypaths  used  in  the  tomography.  Distances  were  restricted  to  1.8  to  15  degrees.  About  50,000 
raypaths  are  shown. 
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3.2  REGULARIZATION  METHOD  AND  PN  ANISOTROPY. 


Tomography  problems  are  always  regularized  in  some  manner.  We  have  chosen  to  use 
Laplacian  damping,  but  there  are  many  viable  methods.  The  damping  parameter  trades-off  the 
image  noise-level  with  image  resolution.  When  damping  is  too  low,  the  resolution  is  good  but 
overwhelmed  by  the  noise;  when  damping  is  too  high  there  is  little  noise  but  the  resolution  width 
is  too  large  to  be  useful.  Between  these  two  damping  extremes  is  an  optimal  level  that  balances 
noise  and  resolution.  When  there  is  only  one  damping  parameter  in  the  problem  the  damping 
level  can  be  set  by  visually  inspecting  the  results  as  is  usually  done  with  seismic  tomography. 
However,  when  multiple  damping  parameters  are  present  visual  inspection  becomes  too 
subjective  to  be  of  use  and  other  criteria  must  be  used.  In  this  paper,  we  are  particularly 
interested  in  the  relative  amount  of  anisotropy  that  is  needed  in  the  solution  and  use  the 
information  criteria  to  resolve  the  issue. 

The  Akaike  Information  Criteria  (AIC)  was  first  introduced  to  determine  the  optimal  model 
order  for  regression  problems  in  statistics  (Akaike,  1973;  McQuarrie  and  Tsai,  1998).  It  also  has 
several  cousins  such  as  the  corrected-AIC  (AICc),  Bayesian  Information  Criteria,  Mallows  Cp- 
criteria,  and  the  Swartz  Information  Criteria  (SIC).  These  criteria  have  been  extended  to 
damped-least  squares  type  problems.  From  a  practical  viewpoint,  they  are  used  to  determine  the 
damping  constants.  We  use  it  to  determine  the  appropriate  levels  of  the  velocity  and  anisotropy 
variations  in  the  seismic  tomography.  From  a  theoretical  viewpoint,  however,  the  information 
criteria  produce  solutions  that  have  optimal  prediction  capability. 

In  this  study  we  present  results  using  the  SIC.  It,  like  the  AIC,  is  implemented  as  a  norm  to  be 
minimized.  For  the  SIC 

SIC  =  log(  SOS  /  n)  +  log(n)  *  k  /  n 

where  SOS  is  the  sum  of  the  squared  residuals  of  the  inversion,  n  is  the  total  number  of  data,  and 
k  is  the  effective  number  of  parameters  of  the  inversion  (McQuarrie  and  Tsai,  1998).  Both  SOS 
and  k  are  functions  of  the  damping  parameters  in  the  inversion.  The  first  term  of  the  SIC  is  the 
logarithm  of  the  variance  of  the  damped  least  squares  inversion.  This  term  decreases  with 
increasing  damping.  The  second  term  is  a  penalty  term  that  decreases  with  increasing  damping. 
It  penalizes  the  solution  for  having  too  many  free  parameters.  Determining  the  minimum  SIC  is 
easily  accomplished  by  running  repeated  damped  least-squares  inversions  for  different  damping 
parameters  and  estimating  the  SIC  for  each  inversion.  The  optimal  damping  parameters  produce 
the  minimum  SIC. 

The  most  difficult  part  of  applying  information  criteria  techniques  to  large  tomographic 
inversions  is  determining  the  number  of  effective  parameters.  For  damped  type  inversions,  the 
trace  of  the  resolution  matrix  gives  an  approximate  estimate  of  the  number  of  effective 
parameters.  However,  the  LSQR  conjugate-gradient  method  used  in  this  study  does  not  produce 
the  resolution  matrix  directly.  We  therefore  use  the  following  approximation  for  the  trace  of  the 
resolution  matrix: 
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k  =  trace  ( (GTG  +  DTD)''  GTG  )  *  trace  (  GTG  (diag  (GTG  +  DTD))'‘ ) 


where  G  is  the  design  matrix  and  D  is  the  damping  matrix.  This  approximation  meets  the  need  of 
being  easily  computed  while  monotonically  varying  between  the  total  number  of  parameters  and 
zero  with  increasing  damping. 

To  choose  the  appropriate  information  norm,  we  tested  several  of  them.  The  AIC,  AICc,  FPE, 
HQ,  and  Cp  norms  were  found  to  substantially  underdamp  the  solutions  while  the  SIC  seemed  to 
slightly  overdamp.  This  is  perhaps  because  the  SIC  was  developed  using  a  Bayesian  criteria  that 
is  more  appropriate  for  large  regression  problems.  A  similar  Bayesian  measure,  the  BIC,  could 
not  be  tested  the  conjugate- gradient  inversion  algorithm  we  use  cannot  compute  the  necessary 
determinants. 

In  our  study,  the  most  important  factor  was  that  the  SIC  gave  an  objective  measure  for  setting  the 
damping  parameters  for  both  the  velocity  and  anisotropy  perturbations.  This  is  particularly 
important  for  determining  the  relative  amount  of  anisotropy  needed  by  the  solution. 

We  used  the  minimum  SIC  to  determine  a  damping  level  of  650  for  the  velocity  parameters 
when  no  anisotropy  was  allowed  in  the  inversion  (Figure  3-2).  At  this  damping  level,  the 
standard  error  is  1 .03  s;  the  rms  level  of  velocity  perturbations  is  0.080  km/s  with  a  maximum 
perturbation  of  0.38km/s.  When  anisotropy  was  included  in  the  inversion  found  the  SIC 
determined  damping  parameters  were  950  and  650  for  the  velocity  and  anisotropy,  respectively 
(Figure  3-3).  The  resulting  inversion  had  an  estimated  standard  error  of  1.07  s;  the  rms  velocity 
variation  is  .047  with  a  maximum  of  0.25  km/s;  the  average  anisotropy  variation  is  .009  km/s  but 
with  a  maximum  of  0.46  km/s.  Bootstrap  estimates  of  all  parameter  errors  found  them  to  be  very 
small. 
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Figure  3-2.  Isotropic  tomography  results  showing  velocity  perturbations  (in  km/s)  relative  to  a  mean  velocity  of  8.0 
km/s.  The  damping  level  determined  by  the  SIC  is  set  at  650. 
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Figure  3-3.  Anisotropic  tomography  results  showing  velocity  and  anisotropy  perturbations  (in  km/s)  relative  to  a 
mean  velocity  of  8.0  km/s.  Much  of  the  travel-time  data  has  been  projected  into  the  anisotropic 
component.  Damping  levels  were  set  using  the  SIC  at  950  and  650  for  the  velocity  and  anisotropy 
components,  respectively. 


The  low  level  of  anisotropy  found  for  most  of  China  makes  it  difficult  to  interpret.  Most  of  the 
anisotropy  is  concentrated  around  Taiwan  and  in  Sichuan  province,  but  results  near  Taiwan  may 
be  heavily  biased  by  poor  locations  and  depth  estimates.  An  anisotropic  checkerboard  test 
(Figure  3-4)  indicates  that  there  is  a  substantial  amount  of  velocity- anisotropy  trade-off  at  the 
chosen  damping  levels.  Specifically,  sharp  variations  in  velocity  map  into  anisotropy.  These 
observations  suggest  that  much  of  the  imaged  anisotropy  is  the  result  of  sharp  variations  in 
velocity.  When  this  same  anisotropic  checkerboard  dataset  was  imaged  without  anisotropy,  the 
results  were  much  better  (Figure  3-5).  Statistically,  the  no-anisotropy  model  has  a  lower  rms 
error  than  the  model  that  includes  anisotropy.  This  again  suggests  that  the  anisotropy,  despite 
the  fact  that  it  fits  the  data  with  fewer  effective  parameters,  is  unnecessary.  The  model  without 
anisotropy  is  our  preferred  model. 
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Figure  3-4.  Checkerboard  test  of  anisotropic  tomography  for  the  SIC  determined  damping  parameters.  The  original 
image  included  both  anisotropy  and  velocity  variations  as  squares  of  alternating  high  and  low  velocity 
with  alternating  east-west  and  north-south  anisotropy.  The  test  used  three-by-three  degree  squares. 
Much  of  the  imaged  anisotropy  has  diagonal  northeast  or  southeast  directions.  This  results  from  the 
velocity  anomaly  trade-off  with  the  anisotropy.  The  figure  shows  velocity  and  anisotropy  perturbations 
(in  km/s)  relative  to  a  mean  velocity  of  8.0  km/s. 
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Figure  3-5.  Checkerboard  test  for  tomography  with  no  anisotropy  but  with  anisotropy  in  the  original  checkerboard 
test  data.  This  test  used  three-by-three  degree  squares.  Other  tests  showed  that  features  as  small  as 
one-by-one  degree  could  be  imaged  in  east-central  china.  The  figure  shows  velocity  perturbations  (in 
km/s)  relative  to  a  mean  velocity  of  8.0  km/s. 


3.3  DISTANCE  DEPENDENCE  OF  THE  RESULTS. 

Lithospheric  structure  may  vary  from  top  to  bottom.  The  rays  used  in  this  study  all  travel 
through  the  lithosphere;  rays  that  dive  into  the  asthenosphere  have  higher  apparent  velocities, 
and  are  not  seen  as  first  arrivals  until  around  20  degrees  distance.  Shorter  raypaths  sample  only 
the  top  of  the  lithosphere  while  more  distant  rays  sample  close  to  the  base  of  the  lithosphere.  To 
investigate  how  important  an  effect  this  is,  the  data  were  divided  into  two  halves  of  nearly  equal 
sizes  and  inverted  tomographically.  Shorter  raypaths,  between  1.8  and  9.5  degrees,  sample  the 
upper-lithosphere.  Longer  raypaths,  between  9.5  and  15  degrees  penetrate  the  top  of  the  upper 
lithosphere  but  are  mostly  in  the  lower  lithosphere.  Comparison  of  the  two  inversions  gives 
insight  into  the  how  the  Chinese  lithosphere  varies  with  depth. 

The  upper-lithosphere  inversion  using  the  shorter  raypaths  is  shown  in  Figure  3-6.  The  mean  Pn 
velocity  for  these  data  was  8.01  km/s.  Western  China  could  not  be  imaged  because  the  station 
spacing  was  too  sparse.  The  image  is  similar  to  that  of  the  full  data  set.  The  long-raypath 
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inversion  (Figure  3-7)  gives  a  similar,  but  different  image  than  the  short-raypath  inversion.  The 
mean  Pn  velocity  was  8.17  km/s  for  these  long  rays  suggesting  that  the  average  mantle  velocity 
varies  by  about  0.2  km/s  from  the  top  to  the  bottom  of  the  lithosphere  at  80-200  km  depth. 
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Figure  3-6.  Inversion  of  short  raypaths  (1.8  to  9.5  degrees  distance).  The  figure  shows  velocity  perturbations  (in 
km/s)  relative  to  a  mean  velocity  of  8.0  km/s. 

The  short-raypath  and  long-raypath  inversions  can  be  interpreted  as  imaging  the  upper  and  lower 
lithosphere.  Similarities  show  features  that  extend  throughout  the  lithosphere  and  differences 
show  features  that  are  predominantly  within  the  upper-  or  lower-lithosphere.  The  major  velocity 
patterns  of  the  two  inversions  are  quite  similar.  In  both  images.  Western  China  has  higher 
velocities  than  in  the  east.  Sichuan  Basin  is  a  prominent  high  velocity  region  that  is  in  both 
images;  the  northeastern  comer  of  Tibet  in  Qinghai  province  is  similar.  This  suggests  that  these 
regions  have  deep  lithospheric  roots  and  may  serve  as  buttresses  to  the  northern  and  eastern 
expansion  of  Tibet.  Differences  between  the  two  images  are  of  interest.  The  rifting  in  Eastern 
China  is  much  more  pronounced  in  the  long-raypath  image,  suggesting  that  these  features  are 
seated  in  the  deeper  lithosphere  and  are  driven  from  the  asthenosphere  below.  Low  velocities 
beneath  Northern  Tibet  are  also  prominent  for  the  long  raypaths,  but  there  was  insufficient  short- 
raypath  data  in  the  area  for  a  comparison.  There  are  several  regions  of  low-velocity  along  the 
eastern  margin  of  Tibet  in  the  short-raypath  images  that  are  not  in  the  long-raypath  image.  These 
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correspond  to  extension  along  the  Panxi  rift  and  Xiaojiang  fault  system.  The  implied 
shallowness  of  these  low-velocity  zones  suggests  that  the  rifting  occurs  primarily  in  the  upper 
lithosphere  is  not  driven  by  asthenospheric  processes.  Rather,  it  is  driven  by  the  oblique 
extension  of  Southeast  Asia  moving  away  from  Tibet  (Royden  and  Burchfiel,  1997;  Chen  et.  al., 
2000). 


Figure  3-7.  Inversion  of  long  raypaths  (9.5  to  15  degrees  distance).  The  figure  shows  velocity  perturbations  (in 
km/s)  relative  to  a  mean  velocity  of  8.0  km/s. 


3.4  TECTONIC  INTERPRETATION. 

The  inversion  yields  an  average  Pn  velocity  of  8.05  km/s  for  all  of  China  consistent  with  the 
worldwide  average  of  8.09  km/s  (Christensen  and  Mooney,  1995).  The  key  results  are  presented 
as  Figures  3-8,  3-9  and  3-10.  Figure  3-8  shows  the  regionally  varying  Pn  velocities 
superimposed  upon  the  topography.  Figures  3-9  and  3-10  are  the  station  delays  and  the  event 
delays.  Both  station  and  event  delays  represent  the  combined  effect  of  crustal  velocity  and 
thickness.  Event  delays,  however,  absorb  any  systematic  errors  in  origin  time  or  depth  in  the 
event  location.  In  most  Pn  tomography  studies  (e.g.  Hearn,  1999],  the  origin-time  and  depth 
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parameters  are  poorly  determined  and  this  creates  such  scatter  in  the  event  delays  that  they  are 
uninterpretable.  This  is  not  the  case  for  the  ABCE  data,  and  it  indicates  that  events  in  China  are 
exceptionally  well  located. 

The  station  delays  can  be  interpreted  as  a  rough  map  of  crustal  thickness  in  China.  For  a 
constant  velocity  crust,  1  second  of  delay  time  represents  approximately  10  km  of  crustal 
thickening.  Late  station  and  event  delays  in  Tibet  then  imply  that  the  crust  there  is  30  to  40  km 
thicker  than  the  crust  in  Eastern  China.  Events  north  of  the  Tarim  Basin  and  in  the  Tien  Shan 
also  show  late  delays  indicating  crust  that  is  10  to  20  km  thicker  than  Eastern  China.  A  thin  crust 
is  also  apparent  southeast  of  Tibet  where  extension  is  occurring. 

The  high  Pn  velocity  regions  are  stable  cratonic  areas  -  Junggar  Basin,  Tarim  Basin,  Taidam 
Basin  and  Sichuan  Basin.  These  cratonic  areas  were  not  greatly  affected  by  tectonic  and 
magmatic  activity  since  its  formation  and  hence  there  is  lower  temperature  in  the  mantle  lid  and 
higher  Pn  velocity.  These  shield  regions  must  also  be  regions  of  strong  lithosphere.  Their  higher 
density  mantle  may  have  aided  in  the  development  of  these  basins.  Pn  velocity  beneath  Tibet 
decreases  from  south  to  north  as  has  been  found  in  previous  Tibet  studies  (Ni  and  Barazangi, 
1983,  1984;  McNamara,  1997).  This  suggests  that  the  subducted  Indian  plate  ends  somewhere 
beneath  central  Tibet.  Eastern  China  has  much  lower  Pn  velocities.  Along  the  Southeastern 
Tibet  margin,  low  Pn  velocities  are  found  suggesting  that  high  temperatures  and  possible  partial 
melt  exists  in  the  mantle  there. 

Variation  in  Pn  velocity  within  the  Tibetan  Plateau  has  been  previously  observed.  Southern 
Tibet  has  a  high  Pn  velocity  and  is  probably  underlain  by  the  under-thrusting  Indian  continental 
lithosphere  (Barazangi  and  Ni,  1983;  Ni  and  Barazangi,  1984;  Beghoul  et  al.,  1993;  Nelson  et  ah, 
1996;  Huang  et  ah,  2000).  Our  data  in  Tibet,  with  limited  ray  coverage,  barely  resolve  the  Pn 
velocity  between  Northern  and  Southern  Tibet.  In  Northern  Tibet,  the  Pn  velocity  is  about  7.9 
km/s  (also  see  McNamara  et  ah,  1997).  The  relatively  low  velocity  in  Northern  Tibet  has  been 
interpreted  as  thermally  weakened  Asian  mantle  lid  that  is  being  squeezed  between  the 
advancing  Indian  lithosphere  to  the  south  and  similarly  rigid  Tsaidam  and  Tarim  lithosphere  to 
the  north  and  west  (Huang  et  ah,  2000).  Our  results  show  that  the  relatively  inactive  Sichuan 
Basin  is  underlain  by  a  deep  high-velocity  lithospheric  root.  It  is  moving  away  from  Tibet  along 
with  the  South  China  Block.  The  Tarim  Basin  and  the  Qinghai  Shan  mountains  also  show  a  deep 
high-velocity  lithosphere,  but  these  are  areas  where  Tibet  is  actively  extending  (Tapponnier, 
2001).  This  high  velocity  mantle  may  have  subducted  beneath  Tibet  from  the  north. 

Beneath  regions  adjacent  to  the  southeastern  edge  of  the  Tibetan  plateau  where  coverage  is  good, 
we  found  low  Pn  velocity  (<7.8  km/s)  regions  in  the  upper  lithosphere  that  do  not  extend  to  the 
base  of  the  lithosphere.  These  regions  correspond  to  magmatic  (extensive  Quaternary  basalt)  and 
tectonic  activity  in  the  Panxi  rift  and  Xiaojiang  fault  system.  The  tendency  of  NNE  alignment  of 
the  Panxi  rift  suggests  an  orientation  controlled  by  the  simple  shearing  motion  in  the  Asian 
lithosphere.  Such  an  interpretation  is  consistent  with  that  the  observation  that  Panxi  rift  and 
Xiaojiang  pull-apart  lie  in  a  trans-extensional  tectonic  environment  produced  by  the  shearing 
between  India  and  Asia  plate,  which  is  currently  taking  place  on  a  transform  boundary  along  the 
right-lateral  Sagaing  fault  (Guzmann-Spazelie  and  Ni,  1993).  These  forces  originate  from  the 
plate  motions  and  takes  place  passively  without  any  asthenospheric  contribution. 
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Eastern  China  is  comprised  of  the  Archean  North  China  block  and  South  China  block.  These 
cratonic  blocks  collided  first  during  the  late  Permian.  After  the  collision,  it  experienced 
widespread  tectonothermal  reactivation  events  during  Late  Mesozoic  and  Cenozoic,  as  indicated 
by  emplacement  of  voluminous  late  Mesozoic  granites  and  extensive  Cenozoic  volcanism 
associated  with  subduction  of  Pacific  and  Philippine  plates  (Yin  and  Nie,  1996).  These 
tectonothermal  events  resulted  in  replacement  of  the  old,  cold,  thick,  and  depleted  lithospheric 
mantle  by  young,  hot,  thin,  and  fertile  lithospheric  mantle  accompanied  by  lithospheric  thinning. 
Widespread  Cenozoic  rifting  and  associated  alkaline  basalt  in  the  North  China  block  is  well 
documented.  GPS  measurements  indicate  that  the  South  China  block  is  moving  southeastward  at 
~1 0  mm/year  while  the  North  China  block  is  moving  eastward  at  ~  5  mm/year  (Chen  et  al.,  2000; 
Michel  et  al.,  2000;  Holt  et  al.,  2000).  The  faster  motion  of  the  South  China  block  produces 
extension  within  the  North  China  block. 


Figure  3-8.  Final  model  of  Pn  velocities  beneath  China.  This  figure  was  constructed  with  raypaths  between  1.8  and 
15  degrees  distance,  nine  arrivals  minimum  per  station  and  event,  and  a  damping  factor  of  650.  The 
figure  shows  velocity  perturbations  (in  km/s)  relative  to  a  mean  velocity  of  8.0  km/s. 
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Figure  3-9.  Final  model  of  Pn  station  delays  beneath  China.  Each  second  of  relative  travel-time  delay  corresponds 
to  approximately  10  km  of  crustal  thickness  change. 

The  Pn  velocity  image  shows  that  low-velocity  structures  are  NNE-elongate  and  correspond  to 
pattern  of  young  volcanism  and  regions  of  large  extension.  Furthermore,  the  low-velocities  are 
more  apparent  in  the  bottom  of  the  lithosphere  suggesting  that  asthenospheric  processes  have 
created  zones  of  partial  melt  above  them.  The  Shanxi  rift,  a  continental  rift  located  on  the 
eastern  edge  of  the  Ordos  Plateau  is  characterized  with  a  very  low  Pn  velocity  that  mimics  the 
en-echelon  pattern  of  the  right  stepping  Shanxi  grabens.  Beneath  the  Tai  Shan,  where  the  Tan 
Lu  fault  crosses,  a  very  pronounced  Pn  velocity  zone  was  imaged.  The  high  elevation  of  this 
mountain  region  must  be  partially  supported  by  the  buoyant  mantle,  a  feature  that  was  not  known 
before  this  study.  The  Hannuoba  alkaline  basalt  that  occurs  north  of  Beijing  and  extends  into 
Inner  Mongolia  is  underlain  by  low  Pn-velocity.  The  overall  low  Pn  velocity  pattern  and 
correlation  with  magmatism  clearly  demonstrate  that  the  lithospheric  mantle  beneath  the  North 
China  block  is  hot  and  partially  melted.  This  interpretation  is  also  consistent  with  the  observed 
inefficient  Sn  propagation  in  this  region  (Rapine  and  Ni,  2002).  The  low-velocities  of  the  North 
China  Block  are  more  apparent  in  the  deeper  lithosphere  (Figure  3-7). 
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The  interior  of  the  North  China  block  consists  of  the  elevated  Ordos  platform  and  the  Yin  Chuan 
rift  located  east  of  the  Ordos  platform.  Pn  velocities  beneath  these  regions  are  normal  and 
indicate  that  Yin  Chuan  rift  is  a  structural  rift  in  an  early  stage  of  rifting.  The  thick  crust  there 
supports  the  Ordos  platform. 

The  eastern  part  of  South  China  block  was  a  magmatic  arc  during  much  of  the  Mesozoic  to 
Cenozoic.  Low  Pn  velocity  regions  correlate  well  with  regions  of  Cenozoic  magmatism.  The 
ray  coverage  near  Taiwan  is  poor  and  thus  imaged  features  are  not  reliable.  A  very  low  velocity 
feature  is  found  north  of  Hainan  Island  in  the  Guangdong  province,  particularly  when  only  the 
long-raypath  data  is  used.  The  close  proximity  of  this  feature  to  the  Red  River  fault  and  South 
China  Sea  suggest  a  possible  relationship  between  the  opening  of  the  South  China  Sea  and  the 
low  Pn  velocity  feature. 
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Figure  3-10.  Final  model  of  Pn  event  delays  beneath  China.  Each  second  of  relative  travel-time  delay  corresponds 
to  approximately  10  km  of  crustal  thickness  change. 
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SECTION  4 


STATION  SPECIFIC  TRAVEL-TIMES 

We  have  prepared  travel-time  curves  for  the  CDSN  stations  in  China  (Figures  4-1).  The  two 
IMS  stations  in  China,  BJT  and  HLA,  were  not  in  existence  when  these  data  were  collected; 
however,  they  replaced  two  nearby  stations,  BJI  and  HLR.  We  determined  Pg,  Pn,  Sg,  and  Sn 
velocities  for  each  of  these  stations  via  straight  line  fits  (Table  4-1).  The  original  phase 
identifications  of  the  ABCE  were  used  with  no  depth  corrections.  The  ABCE  identifies  these 
phases  to  distances  up  to  9  degrees  maximum. 


Table  4-1.  Intercepts  and  apparent  velocities  for  Pg,  Pn,  Sg,  and  Sn  phases  at  CDSN  stations 


Station 

Location 

Pg 

Pn 

Sg 

Sn 

BJI 

Baijiatuan 

-1.1s,  6.14km/s 

4.99s,  7.91  km/s 

-0.56s,  3.66km/s 

10.26s,  4.52km/s 

ENS 

Enshi 

0.8,6.17 

8.6,8.13 

0.5,  3.54 

11.1,4.52 

HLR 

Hailar 

-0.2,  6.06 

6.9,  8.06 

-0.4,  3.52 

10.1,4.51 

KMI 

Kunming 

0.1,6.17 

7.0,  7.94 

0.7,  3.58 

10.9  4.45 

LSA 

Lhasa 

0.3,  6.28 

8.4,  8.02 

0.4,  3.55 

13.5,4.52 

LZH 

Lanzhou 

-1.1,6.07 

8.0,  8.00 

-1.6,3.53 

12.7,  4.55 

MDJ 

Mudanjiang 

-0.6,  6.09 

6.28,8.01 

-1.9,3.52 

9.9,  4.5 

QIZ 

Qiongzhong 

too  little  data 

SSE 

Shanghai 

-0.6,  6.09 

poor  data 

-1.4,3.53 

9.1,4.47 

WMO 

Urumqi 

0.7,  6.29 

8.5,  8.08 

0.7,  3.57 

13.3,  4.55 

numbers  refer  to  intercept  - 

velocity  pairs  in  seconds  and  km/s 

Pn  velocities  surrounding  these  stations  vary  from  7.9  to  8.1  km/s  and  Sn  velocities  are  always 
4.5  km/s.  The  lack  of  variation  in  Pn  and  Sn  velocities  is  due  to  the  large  geographical  averaging 
over  a  circle  of  9  degrees  radius  (254  square-degrees).  Pg  velocities  vary  from  6.0  to  6.3km/s; 
Sg  velocities  vary  between  3.5  to  3.7  km/s.  Higher  Pg  velocities  are  at  Lhasa  and  Urumqi,  both 
regions  of  thicker  crust. 


Intercepts  were  also  computed.  Intercepts  can  be  interpreted  from  the  equation: 
t  =  (2 Z  —  d)J—r - —  where  vc  is  the  crustal  velocity  above  the  refractor,  vn  is  the 

v:  v. 


refractor  velocity,  z  is  the  crustal  thickness,  and  d  is  the  average  event  depth.  The  square-root 
factor  is  approximately  1/10  s/km  for  Pn  and  1/5.6  s/km  for  Sn.  This  was  estimated  using  a  6.2 
km/s  crustal  velocity  and  an  8.0  km/s  mantle  velocity  and  a  conversion  factor  of  0.56  for  the  S- 
wave  velocities.  These  numbers  yield  reasonable  crustal  thicknesses  of  36  to  47  km  if  an 
average  event  depth  of  10  km  is  assumed.  For  Pg  and  Sg,  several  stations  have  negative 
intercepts  for  Pg  and  Sg  as  a  result  of  not  compensating  for  event  depth. 


17 


Time-distance/8.0k/s  Time-distance/8.  Ok/s 


BJI 


ENS 


18 


Time-distance/8.  Ok/s  Time-distance/8.  Ok/s 


HLR 


KMI 


19 


Time-distance/8.  Ok/s  Time-distance/8.  Ok/s 


LSA 


LZH 


20 


Time-distance/8.  Ok/s  Time-distance/8. Ok/s 


MDJ 


QIZ 


21 


SSE 


WMQ 


Figure  4-1.  Reduced  travel-time  plots  for  ten  CDSN  stations.  Stations  BJI  and  HLR  are  next  to  IMS  stations 
BJT  and  HIA.  The  ABCE  routinely  identifies  arrivals  as  Pg,  ,  P,  PP,  pP,  Sg,  Sn,  S,  SS,  sS,  sP, 
PcP,  Sll,  ScS,  PI  1,  ScF,  and  PcS. 
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SECTION  5 


COMPARISON  OF  CHINESE  EVENT  LOCATIONS  WITH  PIDC  LOCATIONS 

We  compared  event  locations  from  the  ABCE  with  those  reported  by  the  Prototype  International 
Data  Center  (PIDC).  These  two  data  sets  overlap  only  for  the  year  1995.  This  is  the  first  year  of 
operation  of  the  PIDC  and  the  last  year  of  data  we  have  for  China.  We  began  by  identifying 
ground  truth  events  for  the  ABCE.  We  used  the  PDIC  criteria  for  defining  ground  truth  events: 
a  GT25  event  was  recorded  by  at  least  50  stations  with  a  maximum  azimuth  gap  of  90°  and  a 
GT10  event  had  at  least  5  stations  with  a  maximum  azimuth  gap  of  180°  within  3-degree  distance 
(Yang  and  Romney,  CMR  Technical  Report  99-15,  1999).  These  ground  truth  events  were 
compared  to  locations  from  the  PIDC  Relocated  Event  Bulletin  (REB).  For  the  GT10  events,  23 
of  91  events  were  reported  by  the  PIDC  (Figure  5-1).  For  the  GT25  events,  86  of  89  events  were 
reported  by  the  PIDC  (Figure  5-2).  Events  that  were  not  reported  by  the  PIDC  were  generally 
less  than  magnitude  4.  On  average,  mislocations  are  less  than  22  km  distance,  but  several  events 
are  mislocated  by  over  33  km.  The  largest  location  difference  was  by  1 18  km  for  a  magnitude 
3.6  event  on  October  30,  1995  (8:26)  in  Yunnan  province.  For  this  event,  the  PIDC  location 
used  5  stations  with  a  210°  azimuth  gap  and  lists  standard  errors  of  between  50  to  60  km.  The 
PIDC  depth  was  constrained  to  0  km.  The  National  Earthquake  Information  Center  (NEIC) 
reported  the  event  101  km  away  from  the  ABCE  location  and  the  International  Seismological 
Centre  (ISC)  reports  the  event  17  km  away.  Curiously,  the  location  given  to  the  ISC  by  Beijing 
was  off  by  14  km  from  the  ABCE  location  that  we  obtained. 
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Figure  5-1. 


Relative  locations  from  23  GT10  events  reported  by  both  the  PIDC  and  the  Annual  Bulletin  of  Chinese 
Earthquakes  1995.  Axis  scales  are  in  degrees. 
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Figure  5-2.  Relative  locations  from  86  GT25  events  reported  by  both  the  PIDC  and  the  Annual  Bulletin  of  Chinese 
Earthquakes  1995.  Axis  scales  are  in  degrees. 


There  were  a  limited  number  of  ground  truth  events,  so  the  entire  1995  ABCE  catalog  was 
searched  for  all  events  located  by  both  the  ABCE  and  PDIC.  Histograms  of  earthquake  size  in 
China  for  1995  are  plotted  for  both  the  ABCE  and  the  PIDC  in  Figure  5-3.  These  show  that  the 
ABCE  locates  slightly  more  events  in  all  magnitude  ranges  than  the  PIDC.  This  discrepancy  is 
from  differences  in  the  magnitude  scales  used  by  ABCE  (local  magnitude)  and  PIDC  (body- 
wave  magnitude).  The  plot  deviates  from  linear  for  magnitudes  less  than  4.0  to  4.5  for  both 
networks  suggesting  that  smaller  events  are  often  missed  by  both  catalogs. 

We  found  195  events  that  were  located  by  both  the  ABCE  and  PIDC  during  1995. 
Mislocation  vectors  for  these  are  shown  in  Figures  5-4  and  5-5.  Many  of  the  events  are  in 
Taiwan,  where  the  ABCE  locations  may  not  be  as  accurate.  Differences  between  the  ABCE 
epicentral  locations  and  the  equivalent  PIDC  locations  average  0.95  degrees  with  a  median  of 
0.55  degrees  (Figure  5-6).  The  sample  is  biased  by  a  few  events  with  large  location  differences 
of  up  to  four  degrees.  Events  with  extreme  location  differences  of  over  three  degrees  were  found 
to  have  limited  station  coverage,  with  five  or  fewer  stations  recording,  and  azimuthal  gaps  larger 
than  1 80  degrees.  For  the  events  in  our  sample,  the  estimated  event  location  errors  listed  in  the 
PIDC  catalog  averaged  ±0.76  degrees.  Considering  the  outliers  in  Figure  5-6,  the  PIDC  error 
estimates  are  reliable.  Poorly  located  events  can  be  identified  on  the  basis  of  their  estimated 
error  ellipses  and  the  station  coverage. 

The  CTBT  requires  that  events  be  located  within  a  1000  km-squared  region  at  the  2- 
sigma  (95%)  confidence  level.  The  results  find  that  this  was  not  possible  for  the  PIDC  for  most 
events  during  1 995  because  of  inadequate  station  coverage.  This  is  also  reflected  in  the  large 
error  estimates  in  the  PIDC  catalog.  The  best  method  to  improve  event  location  in  China  is  to 
improve  station  coverage. 
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curve  of  PIDC  and  ABCE 


Figure  5-3.  Histogram  of  event  counts  in  different  magnitude  ranges  for  the  Annual  Bulletin  of  Chinese 
Earthquakes  (ABCE)  and  the  Preliminary  International  Data  Center  (PIDC)  for  1995.  The  lower 
completeness  threshold  is  between  magnitude  4.0  and  4.5;  some  of  the  differences  are  due  to 
differences  in  magnitude  determination  between  the  two  catalogs. 
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Figure  5-4.  Map  showing  location  differences  between  the  ABCE  and  PIDC  catalogs  for  China  for  1995.  Arrows 
begin  at  the  ABCE  location  and  end  at  the  PIDC  location. 
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Figure  5-5.  Composite  plot  of  195  location  difference  vectors  between  ABCE  and  PIDC  for  China  1995.  Vectors 
point  from  the  ABCE  location  to  the  PIDC  location.  Axis  scales  in  degrees  distance. 
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Figure  5-6.  Histogram  of  location  differences  between  the  ABCE  and  PIDC  locations  for  China,  1995.  The 
average  location  difference  is  0.95  degrees;  however,  this  is  biased  by  a  few  poorly  located  events.  The 
median  location  difference  is  0.55  degrees  and  the  most  probable  location  difference  is  0.1  degrees. 
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SECTION  6 


PROPAGATION  CHARACTERISTICS  OF  SN  AND  LG  IN  NORTHEASTERN  CHINA 

AND  MONGOLIA 

We  present  propagation  characteristics  of  the  short-period  seismic  phases  Sn  and  Lg  in 
Northeastern  China  and  Mongolia.  Richard  Rapine  and  James  Ni  have  submitted  this  section  of 
the  report  to  the  Bulletin  of  the  Seismological  Society  of  America.  Three-component  broadband 
and  short-period  seismograms  from  the  Global  Seismic  Network  station  ULN  and  the  Chinese 
Digital  Seismic  Network  stations  BJT,  HIA,  and  MDJ  were  used  in  this  study.  The  digital 
waveforms  were  bandpass  filtered  (1. 0-5.0  Hz)  and  were  used  to  rank  Sn  and  Lg  amplitudes 
relative  to  the  Pn  and  Pg  phases.  The  short-period  Sn  waveforms  exhibit  great  variability. 
Efficient  Sn  propagation  is  observed  in  Southern  Siberia  and  along  the  Daxinganling  Mountains. 
Inefficient  Sn  transmission  is  found  near  Lake  Baikal  and  most  of  Northeastern  China.  Sn  does 
not  propagate  across  Mongolia  or  the  North  China  Plain.  In  Northeast  China,  zones  of 
attenuation  exist  within  regions  of  extension  and  high  heat  flow.  In  Mongolia,  severe  Sn 
attenuation  suggests  the  presence  of  a  partially  melted  uppermost  mantle.  The  zones  of 
inefficient  Sn  propagation  coincide  with  zones  of  low  Pn  velocities  (7.9  km/s)  determined  by  Wu 
(1998).  In  contrast,  the  Lg  seismic  phase  has  large  amplitudes  and  propagates  efficiently  for  all 
distances  in  Northeastern  China  and  Mongolia. 

Heterogeneity  and  anelasticity  of  the  crust  and  upper  mantle  cause  amplitude  reduction  of 
regional  seismic  phases.  Measurements  of  this  loss  in  amplitude  can  be  used  to  infer  the 
physical  properties  and  structure  of  a  particular  region  on  earth.  The  propagation  of  short-period 
regional  Sn  and  Lg  waves  in  Northeastern  China,  Mongolia,  and  Southern  Siberia  is 
characterized  here  and  used  to  investigate  how  tectonic  and  geologic  structures  affect  these 
regional  phases.  This  study  utilizes  the  digital  broadband  and  short-period  data  from  the  Global 
Seismic  Network  (GSN)  station  ULN  (Ulaanbaatar,  Mongolia)  and  the  Chinese  Digital  Seismic 
Network  (CDSN)  stations  BJT  (Beijing,  China),  HIA  (Hailar,  China),  and  MDJ  (Mudanjiang, 
China).  Propagation  efficiencies  of  Sn  and  Lg  phases  were  qualitatively  analyzed  by  ranking 
their  amplitudes  relative  to  the  Pn  and  Pg  phases.  The  data  were  then  used  to  determine  the  areas 
of  Sn  and  Lg  attenuation  in  the  region  of  Northeastern  China  and  Mongolia. 

The  North  China  Plain  and  the  Songliao  Basin  of  Northeastern  China  are  large  Cenozoic 
extensional  structures  (Ren  et  al.,  1987).  The  locations  of  these  basins  coincide  with  areas  in 
which  the  lithosphere  is  thinned  and  high  heat  flow  is  present.  These  basins  represent  the  effects 
of  Mesozoic-Cenozoic  (Songliao  Basin)  and  Cenozoic  crustal  extension  (Western  North  China 
Plain)  proceeded  by  Andean-type  subduction,  related  deformation  and  magmatism.  The 
magmatic  arc  during  the  Mesozoic  subduction  episode  modified  the  Daxinganling  Mountains.  It 
has  been  suggested  that  the  cause  of  extension  was  the  upward  advection  of  the  asthenosphere, 
which  was  related  to  the  foundering  of  the  Kula-Pacific  plates  beneath  the  Eastern  Asia 
continental  margin  since  the  late  Mesozoic  time  (Ren  et  al.,  1987). 

Variations  in  the  structure  and  properties  of  the  crust  and  uppermost  mantle  greatly  affect  short- 
period  seismic  waves  recorded  at  regional  distances  (200-2000  km).  Lg  waves,  which  are 
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composed  of  shear  waves  that  have  been  trapped  in  the  crust,  have  characteristics  that  are 
diagnostic  of  the  crustal  structure.  In  contrast,  Sn  waves  travel  in  the  mantle  lid  and  are 
indicative  of  the  physical  properties  there.  Changing  of  signal-to-noise  amplitude  ratios  of  these 
regional  phases  can  provide  insight  into  where  these  variations  occur.  Since  regional  seismic 
phases  are  used  for  discriminating  nuclear  explosions  from  earthquakes  (e.g.,  Taylor  et  al., 
1989),  it  is  important  to  know  where  these  phases  are  attenuated.  Major  structural  features  such 
as  sedimentary  basins,  rapidly  changing  surface  and  Moho  topography  will  all  influence  the 
character  of  regional  seismograms  (e.g.,  Gibson  and  Campillo,  1994).  It  is  clearly  important  to 
map  where  these  features  affect  the  propagation  of  regional  seismic  phases.  This  kind  of 
regional  data  is  an  integral  part  of  the  effort  toward  verification  of  a  Comprehensive  Test  Ban 
Treaty  (CTBT). 

Figure  6-1  shows  the  seismic  station  and  event  locations  used  in  this  study.  These  stations  are 
part  of  the  Chinese  Digital  Seismic  Network  (CDSN)  and  Global  Seismic  Network  (GSN). 
Table  6-1  summarizes  the  station  information  and  includes  the  years  of  data  we  analyzed  from 
each  station.  Short-period  data  were  obtained  from  stations  BJT,  HIA,  and  MDJ  in  China. 
Broadband  data  were  obtained  from  ULN  in  Mongolia.  The  broadband  data  were  recorded  on  a 
three-component  Streckeisen  STS-1  instrument.  The  STS-1  records  data  at  20  samples/s  and  has 
a  flat  amplitude  response  from  .01  to  10  Hz.  Short-period  data  were  recorded  by  Model  DJ-1 
seismometers  at  40  samples/s.  The  seismic  events  requested  from  the  Incorporated  Research 
Institutions  for  Seismology  Data  Management  Center  (IRIS  DMC)  had  body-wave  magnitudes 
greater  than  4.3,  epicentral  distances  between  2  and  15  degrees,  the  complete  range  of 
backazimuths,  and  depths  less  than  50  km.  Event  origin  times  and  locations  were  taken  from  the 
Preliminary  Determination  of  Epicenters  (PDE)  catalogs  distributed  by  the  U.S.  Geological 
Survey  (USGS).  All  data  with  low  signal-to-noise  ratios  were  removed  from  the  data  set  used 
for  this  analysis. 
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Figure  6-1.  A  station  and  event  location  map  for  the  region  of  study.  Triangles  represent  CDSN  and  GSN  station 
locations.  Circles  represent  event  locations  and  are  proportional  to  the  event  magnitudes. 


Short-period  and  broadband  data  were  used  to  study  the  regional  wave  propagation  of  Sn  and  Lg 
phases  in  Northeast  China  and  Mongolia.  The  high-frequency  propagation  efficiencies  of  these 
phases  were  used  to  learn  more  about  the  structure  and  composition  of  the  crust  and  upper 
mantle  in  the  area  of  interest.  This  approach  has  been  used  many  times  before  for  different 
geologic  regions  (e.g.,  Kadinsky-Cade  et  al.,  1981;  Rapine  et  al.,  1997;  Rodgers  et  al.,  1997).  A 
bandpass  filter  was  used  to  filter  the  data  between  1. 0-5.0  Hz. 


Table  6-1.  Summary  of  the  stations,  recording  period  of  data  and  type  of  data. 


Station  Name 

Lat(°N) 

Lon(°E) 

Data  Period 

Type* 

BJT 

40.10 

116.17 

1994-1997 

SP 

HIA 

49.27 

119.74 

1988-1996 

SP 

MDJ 

44.62 

129.59 

1989-1996 

SP 

ULN 

47.87 

107.05 

1994-1997 

B 

*  SP  =  short-period,  B  =  broadband 


Sn  waves  predominantly  travel  within  the  uppermost  mantle  and  since  they  are  shear  waves,  they 
are  more  sensitive  to  attenuation  than  P  waves.  At  regional  distances,  Sn  arrives  with  a  group 
velocity  of  4. 4-4.7  km/s.  Sn  propagation  efficiencies  have  been  used  to  infer  the  physical 
properties  of  the  mantle  lid.  For  example,  Sn  generally  does  not  propagate  across  regions  of  high 
heat  flow  such  as  the  Basin  and  Range  Province  and  the  central/northem  Tibetan  Plateau  (e.g., 
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Ni  and  Barazangi,  1983;  Beghoul  et  al.,  1993;  McNamara  et  al.,  1995).  Poor  Sn  propagation 
may  be  attributed  to  high  shear  wave  attenuation  associated  with  a  partially  melted  mantle  lid 
(Molnar  and  Oliver,  1969).  Efficient  Sn  propagation  is  observed  beneath  stable  shield  regions 
(e.g.,  Brune  and  Dorman,  1963;  Bath,  1966). 

Lg  waves  are  a  complex,  high-frequency  wavetrain  with  large  amplitudes  at  regional  distances 
and  group  velocities  from  3.2  to  3.7  km/s.  Lg  has  been  modeled  as  a  composite  of  higher-mode 
Love  and  Rayleigh  waves.  In  terms  of  ray  theory,  Lg  can  be  considered  as  a  sequence  of 
multiply  reflected,  postcritical  S  waves  trapped  in  the  continental  crust  (Bouchon,  1982).  Lg 
propagates  efficiently  in  continental  shields  but  is  not  seen  across  ocean  basins  due  to  the  thin 
oceanic  crust  (Press  and  Ewing,  1952).  The  amplitude  of  Lg  may  be  weakened  by  variations  in 
crustal  structure  along  the  propagation  path  or  by  a  low-Q  or  partially  melted  crust  (Ruzaikin  et 
al.,  1977;  Campillo,  1987;  Husebye  and  Ruud,  1996). 

The  high  signal-to-noise  data  were  used  to  manually  pick  the  regional  phase  arrivals  of  Sn  and 
Lg.  We  picked  the  arrival  times  of  Sn  and  Lg  to  coincide  with  group  velocities  of  4.5  km/s  and 
3.5  km/s,  respectively.  Previous  studies  analyzed  short-period  seismograms  to  qualitatively 
characterize  broad  tectonic  regions  in  terms  of  efficient  or  inefficient  Pn,  Sn,  and  Lg  propagation 
(e.g.,  Ruzaikin  et  al.,  1977;  Kadinsky-Cade  et  al.,  1981;  Ni  and  Barazangi,  1983;  Rodgers  et  al., 
1997).  We  used  a  similar  approach  to  rank  the  efficiencies  of  the  observed  regional  phases.  We 
classified  efficient  Sn  and  Lg  phases  as  those  phases  with  signal-to-noise  (S/N)  ratios  greater 
than  or  equal  to  two.  Sn  and  Lg  phases  with  S/N  ratios  between  one  and  two  were  classified  as 
inefficient.  When  no  Sn  or  Lg  phase  can  be  seen  on  a  seismogram,  we  classified  that  phase  as 
not  propagating  for  that  station-event  path. 

6.1  RESULTS. 

Representative  short-period  seismograms  recorded  by  CDSN  stations  are  provided  in  Figures  6-2 
and  6-3.  Efficient  Sn  propagation  is  observed  in  Southern  Siberia  and  the  Daxinganling  (Figure 
6-4).  Apparently,  the  lithosphere  beneath  this  mountain  is  sufficiently  cooled  and  does  not 
contain  partial  melt.  Poor  Sn  transmission  is  found  in  Lake  Baikal  and  most  of  Northeastern 
China  where  late  Mesozoic  to  Cenozoic  extension  has  occurred  (Figure  6-5).  Sn  does  not 
propagate  across  Mongolia  and  most  of  the  North  China  Plain,  where  Late  Tertiary  to 
Quaternary  volcanism  exists  (Figure  6-6).  Rapine  et  al.  (1997)  also  observed  no  Sn  propagation 
across  Mongolia  from  earthquakes  that  occurred  in  Siberia  and  recorded  at  stations  in  China. 
The  strong  Sn  attenuation  beneath  the  North  China  Plain  and  Mongolia  suggests  that  there  is 
partial  melt  in  the  uppermost  mantle  beneath  those  regions. 

High-frequency  Lg  waves  propagate  efficiently  in  Southern  Siberia,  Northeastern  China,  and 
Mongolia  (Figure  6-7).  In  this  region  of  study,  there  were  no  observations  of  Lg  blockage  or 
attenuation.  Lg  amplitudes  remained  large  for  even  the  longest  distances  studied.  Figures  6-2 
and  6-3  show  the  large  Lg  amplitudes  relative  to  other  phases  and  the  noise.  Structural 
boundaries  in  the  crust  in  this  region  do  not  have  a  large  effect  on  Lg  propagation  as  they  do  in 
such  regions  as  the  Tibetan  Plateau. 
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BJT  SHN  JAN  5  (005),  1 997  Delta  =  9.46°  BAZ  =  246.93° 


HIA  SHE  JAN  13  (013),  1 991  Delta  =13.16°  BAZ  =  233.80° 


Figure  6-2.  Representative  seismograms  from  a  short-period  CDSN  station.  A  time  axis  is  provided  at  the  top  of 
the  Figure,  zero  corresponds  to  the  Pn  arrival  time.  Vertical  bars  represent  the  expected  arrival  time  of 
Sn  or  Lg.  All  seismograms  were  bandpass  filtered  between  1. 0-5.0  Hz. 
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HIA  SHE  JUL  18(199),  1990  Delta  -  7.99°  BAZ  =  323.96° 


HIA  SHE  OCT 25  (298),  1989  Delta=8.29°  BAZ  =  356.51° 
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HIA  SHE  FEB  25  (056),  1983  Delta  =  7.20°  BAZ  =  1 63.56° 


Figure  6-3.  Representative  seismograms  from  short-period  and  broadband  stations  of  the  CDSN  and  GSN,  same  as 
Figure  2. 
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Figure  6-4.  A  map  of  efficient  Sn  wave  propagation  in  China,  Mongolia  and  Southern  Siberia.  Triangles  represent 
station  locations,  circles  represent  event  locations  and  the  numbers  refer  to  seismograms  from  Figures  2 
and  3. 


90*  100*  110*  120*  130*  140* 


Figure  6-5.  A  map  of  inefficient  Sn  wave  propagation  paths.  Triangles  represent  station  locations,  circles  represent 
event  locations  and  the  numbers  refer  to  seismograms  from  Figures  2  and  3. 
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90*  100’  110*  120*  130*  140* 


Figure  6-6.  A  map  shows  the  station-event  paths  where  no  Sn  propagates.  Triangles  represent  station  locations, 
circles  represent  event  locations  and  the  numbers  refer  to  seismograms  from  Figures  2  and  3. 
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Figure  6-7.  A  map  of  efficient  Lg  wave  propagation  in  China,  Mongolia  and  Southern  Siberia.  Triangles  represent 
station  locations,  circles  represent  event  locations  and  the  numbers  refer  to  seismograms  from  Figures  2 
and  3. 
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Figure  6-8.  A  summary  map  of  efficient,  inefficient  and  no  Sn  wave  propagation  in  Northeast  China  and 
surrounding  regions. 


6.2  DISCUSSION. 

The  propagation  efficiencies  of  Sn  and  Lg  waves  have  been  mapped  for  Northeast  China  and 
Mongolia.  A  summary  map  of  the  gross  characteristics  of  short-period  Sn  regional  waves  has 
been  provided  in  Figure  6-8.  This  map  provides  some  insight  into  the  uppermost  mantle 
properties  of  the  region  and  explicitly  shows  the  regions  of  attenuation  we  observed.  Lg 
propagates  efficiently  and  shows  large  amplitudes  for  all  paths  studied  in  the  region.  The 
Cenozoic  tectonics  of  Northeast  China  has  had  an  influence  on  short-period  Sn  wave 
propagation.  A  P  wave  study  of  the  upper  mantle  beneath  Northern  China  and  Mongolia 
measured  low  P  wave  velocities  in  the  uppermost  mantle  (7.9  km/s)  in  most  of  the  region  of  our 
study  (Wu  et  al.,  1998).  This  Sn  study  provides  evidence  that  the  uppermost  mantle  may  be 
partially  melted  in  parts  of  Mongolia  and  Northeast  China.  A  buoyant  mantle  beneath  the 
Mongolian  Plateau  would  support  the  high  topography  observed  there.  Any  partially  melted 
mantle  is  probably  a  result  of  either  convective  removal  of  part  of  the  mantle  lid  or  prior 
subduction  (Paleozoic)  of  water-rich  materials  into  the  upper  mantle  leading  to  a  present 
weakening  and  melting  due  to  strain  heating.  This  strain  heating  is  probably  a  result  of  the 
Indian-Eurasian  collision  whose  effects  have  been  felt  all  the  way  up  to  Siberia  (e.g.,  Tapponnier 
and  Molnar,  1977). 
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SECTION  7 


LATERAL  VARIATION  OF  LG  ATTENUATION  IN  TIBETAN  PLATEAU 

7.1  INTRODUCTION. 

Lg  is  a  crustal  guided  wave  propagating  efficiently  for  several  thousands  of  kilometers  in  stable 
continental  shields  with  the  group  velocities  between  3.7  and  3.0  km/s.  Modeling  of  Lg  has 
shown  that  it  is  the  sum  of  higher  mode  surface  waves  or  superposition  of  multiply  reflected 
shear  waves  within  the  crust  (Bouchon,  1982;  Kennett,  1986).  The  characteristic  of  Lg  phase  is 
sensitive  to  changes  within  the  crustal  wave-guide  such  as  sedimentary  basins  and  large  lateral 
variations  in  depth  to  Moho  and  in  crustal  seismic  velocity  structure  (Ruzaikin  et  al.,  1977; 
Campillo,  1987).  Press  and  Ewing  (1952)  found  that  Lg  does  not  propagate  within  thin  oceanic 
crust,  while  the  presence  of  Lg  phase  is  often  associated  with  path  traveling  within  continent 
(Oliver  et  al.,  1955).  Lg  attenuation  is  caused  by  two  main  sources:  scattering  due  to 
heterogeneities  and  intrinsic  absorption.  Many  attempts  have  been  made  to  separate  these 
contributions  (e.g.  Mayeda  et.  al.,  1992).  Their  results  indicate  that  the  proportion  of  each 
contribution  might  be  different  for  different  regions. 

Tibet  is  the  one  of  the  most  striking  example  of  continent-continent  collision.  The  collision 
began  about  50  Ma  causing  the  Tibetan  Plateau  rising  to  its  5  km  elevation.  Although  Lg  is  one 
of  the  most  dominant  phases  observed  on  regional  seismograms  for  paths  traveling  within 
continental  shield,  studies  in  Tibet  shows  that  the  maximum  distance  to  observe  Lg  wave  is  about 
600  km  within  the  Plateau  (McNamara,  1996).  Beyond  that  distance,  Lg  energy  is  largely 
attenuated.  There  have  been  a  number  of  studies  to  find  Lg  attenuation  in  Tibet.  Ni  and 
Barazangi  (1983)  qualitatively  mapped  the  Lg  attenuation  where  they  found  inefficient  Lg  in  the 
plateau.  Rapine  et  al.  (1996)  found  that  Lg  is  generated  in  Tibet,  however,  paths  to  LSA  (Figure 
7-1)  exhibits  large  attenuation.  They  also  suggested  that  the  plateau  boundaries  causes  scattering 
of  Lg  energy.  Although,  high  attenuation  in  Tibet  is  well  known,  none  of  the  methods  could  find 
Lg  Q0  values  much  less  than  200  for  Tibetan  Plateau  (Reese  et  al.,  1999).  Xie  (2002)  used  two- 
station  method  allowing  backazimuthal  gap  up  to  36°  and  found  Lg  Qo  =126  for  Eastern  Tibet. 
He  averaged  all  two-station  measurements  to  obtain  an  average  Lg  Q.  Because  high  frequency 
Lg  is  rapidly  attenuated  in  areas  of  very  low  Q,  it  is  especially  difficult  to  quantify  areas  of  low 
Q  such  as  those  in  Southernmost  Tibet  because  the  lack  to  stations  in  Tibet.  In  this  study,  we  use 
data  from  several  experiments  conducted  within  Tibet  and  applied  two  different  techniques  to 
determine  the  lateral  variation  of  Lg  Q  within  Tibet  (Gok  et  al.,  2004).  We  also  mapped  the 
presence  of  lateral  change  along  the  profiles  of  INDEPTH  II  and  INDEPTH  III. 

7.2  DATA  AND  METHOD. 

During  the  period  of  1991-1999  there  were  three  PASSCAL  temporary  broad-band  seismic 
experiments  to  explore  the  crustal  and  upper  mantle  structure  of  Tibet  (Figure  7-1).  The  first  one 
was  the  1991-1992  IRIS-PASSCAL  experiment,  a  one- year  deployment  in  the  eastern  part  of 
Tibet  with  1 1  broad-band  stations.  The  second  one  was  INDEPTH  II  extending  from  the  High 
Himalayas  to  about  150  km  north  of  the  IYS  (Indus- Yarlung  Suture).  Nine  of  the  fifteen  stations 
equipped  with  broad-band  CMG-3T  seismometers.  The  remaining  six  stations  were  equipped 
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with  1  Hz  short  period  seismometers.  The  third  experiment,  ENDEPTH  III,  operated  for  one  year 
(July,  1998-June,  1999)  across  the  Bangong-Nujiang  Suture  (BNS)  with  more  stations  (37 
broadband,  10  semi-broadband  and  15  short-period  stations).  INDEPTH  II  and  INDEPTH  III 
provided  a  unique  opportunity  to  measure  the  lateral  variation  of  Lg  attenuation  in  Tibet. 
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Figure  7-1.  Station  distribution  of  seismic  experiments  during  period  of  1991-1999. 


37 


Station:BB05  Starting  time  =  106.6sec.  Dist  =4  3  7 .0  l<ni 
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Figure  7-2.  Lg  wave  package  and  noise  part  spectra  of  event  1994  236. 


We  analyzed  ten  earthquakes  from  INDEPTH  II,  eleven  from  1991-1992  PASSCAL  experiment 
and  six  from  INDEPTH  III.  All  selected  earthquakes  have  crustal  depths.  We  used  only  vertical 
component  seismograms  within  distances  of  2°  to  6°.  The  earthquake  magnitudes  range  between 
4.0  <  Mb  <  5.5.  The  Lg  group  velocity  window  is  set  between  3. 7-3.0  km/s  depending  on  the 
wave  package  (Figure  7-2).  In  this  study,  we  used  two  different  frequency  domain  techniques  to 
estimate  Lg  attenuation  in  Tibet.  Both  methods  have  been  extensively  tested  on  earthquakes  and 
explosions  (Xie  1993,  1998). 


7.3  STOCHASTIC-SPECTRAL  INVERSION  OF  LG  Q0. 


Nuttli  (1973)  first  used  Lg  amplitudes  in  the  time  domain  to  measure  path  attenuation  and  the 
size  of  the  seismic  source.  In  frequency  domain,  the  spectral  inversion  method  of  Xie  (1993) 
was  developed  to  simultaneously  invert  for  source  and  path  variable  Q  unlike  the  other  inversion 
methods  which  were  based  on  either  an  assumption  of  constant  Q  or  separately  measured  source 
and  attenuation  terms.  Lg  displacement  amplitude  spectra  can  be  expressed  as  a  product  of  terms 
in  frequency  domain  at  i  station, 
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(7.1) 


using  the  frequency  dependency  of  Lg, 


Q(f)  =  QoP 


(7.2) 
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where  S(f)  is  the  source  spectrum  assuming  Lg  is  not  effected  by  the  radiation  pattern  and  can  be 
written  as, 


S(f)  = 


Mp  1 
Anpv]  l  +  f2/f2 


(7.3) 


where  p  and  v  are  the  averaged  crustal  density  and  shear-wave  velocity,  respectively  (Street  et 
al.,  1975).  GO)  is  the  geometrical  spreading  term,  y[\Aj  ,  where  A0=100  km.  X(f)  is  site 

response  and  R(f)  is  randomness  due  to  propagation  such  as  random  effects  of  multiple  ray  Lg 
arrivals  and/or  3D  heterogeneity.  Using  the  above  relation  we  end  up  with  the  model  vector 

mT  =(M0,fc,  Qov  .,MN,fN,Q0N,TJN)r  (7.4) 

The  radiation  pattern,  systematic  deviation  in  velocity  structure,  station  site  response  and 
systematic  deviation  from  w'2  model  is  ignored.  Equation  (7.1)  is  linearized  by  taking  the 
logarithm  of  both  sides.  For  a  single  event  recorded  by  N  stations,  unknowns  are  estimated  by 
grid  searching  over  the  possible  Mo  and  fc  using  the  following  relationship, 
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i  =  1,2,..., N  (7.5) 


Qoi  and  Oj  are  estimated  over  this  loop,  and  the  complete  model  vector  is  used  to  estimate  a 
residual  square, 
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The  residual  square  is  estimated  for  the  all  possible  Mo  and  fc  values.  The  model  vector  that 
minimizes  the  solution  is  the  best  solution  (Xie,  1998). 

7.4  TWO-STATION  INVERSE  METHOD. 


The  two-station  spectral  ratio  method  is  a  very  robust  technique  of  measuring  attenuation  since  it 
properly  eliminates  source  effects  and  results  in  a  less  biased  estimate  of  Lg  Qo.  In  this  method, 
the  pair  of  stations  is  selected  so  that  they  lie  on  the  same  great  circle  path  as  from  the  source. 
The  Lg  Q(f)  in  the  study  area  is  assumed  as  being  laterally  homogeneous.  When  an  event  is 
recorded  by  two  stations,  j  1  and  j2,  along  the  same  great  circle  path,  we  denote  the  Lg  amplitude 
spectra  recorded  at  the  station  j  as  Aj  (j  =jl  j2),  and  define  a  scaled  logarithmic  spectral  ratio 
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where  Vg  is  the  Lg  group  velocity,  R(jl)  and  R(j2)  are  the  epicentral  distances.  For  each  of  the 
station  pairs  (j  1,  j2),  the  number  of  events  that  are  recorded  along  the  same  great  circles  as  I  (j  1 , 
j2).  The  mathematical  expectation  was  given  by  (Xie  and  Mitchell  1990;  Xie,  2002) 


E(D,)  = 


(7.8) 


(7.9) 


7.5  LATERAL  VARIATION  OF  LG  ATTENUATION  . 


The  Algebraic  Reconstruction  Technique  (ART)  is  a  way  of  imaging  lateral  changes  of  velocities 
or  amplitudes.  Xie  and  Mitchell  (1990)  applied  an  imaging  algorithm  that  is  originally  for  the 
regionalization  of  Lg  coda  Q  estimates.  We  adopted  this  technique  for  use  here.  Since  we  aim  to 
resolve  only  the  change  along  the  array  we  used  only  the  ID  imaging. 


The  ART  technique  can  be  briefly  summarized  as  follows, 
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Figure  7-3.  We  assumed  a  1-D  profile  for  Q  inversion.  The  dashed  line  is  the  great  circle  connecting  stations  1 
and  2.  Inter-station  Q  is  measured  for  the  dashed  line  path.  All  possible  two-station  pairs  are  used  for 
the  inversion. 


Smn  is  the  area  where  the  nth  path  overlaps  with  the  mth  grid,  Sn  is  the  total  cells  for  each  path, 
QJ,,  is  a  starting  model  for  m  being  m=l,2,3,...,Ng  (Figure  7-3).  For  each  of  the  m  grids 

(m=l,2,3,...,Ng),  a  new  estimate  of  Qm,  denoted  by  Q1*1 ,  is  made  by  back  projecting  A'(I  into  the 
inverse  of  the  Q'm  ’s 
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To  stabilize  the  convergence  and  smooth  the  model  Suetsugu  and  Nakanishi  (1985)  proposed  a 
nine  point  spatial  smoothing  procedure  (Xie  and  Mitchell,  1990,  Equation  14).  We  used  less 
smoothing  points  then  the  original  method  since  we  are  only  resolving  1 D  variation. 

7.6  ERROR  ANALYSIS. 

The  error  analysis  can  be  approximated  using  standard  deviation  of  each  single  trace  Q 
measurements.  It  is  assumed  that  this  ctq  gives  a  measure  of  the  absolute  value  of  statistical  error 
(Xie  and  Mitchell,  1990).  The  error  values  are  empirically  constructed  using  the  mean  of  <7q. 
This  mean  value  is  used  to  generate  Gaussian  random  noise  series  whose  number  is  20xNg.  The 
number  of  generated  noise  series  are  then  added  to  Qn, 

e;= q,+sq,  (7.i2) 

These  Qn  values  are  inverted  20  times  and  compared  with  original  Qm  image.  The  difference 
gives  the  amount  of  error  in  coda  Q  estimate. 

Back  projection  tomography  does  not  yield  a  resolution  matrix.  Xie  and  Mitchell  (1990)  used  a 
technique  that  is  based  on  an  approximation  to  the  resolution  kernel.  By  reconstructing  the 
chosen  block  with  a  unit  value,  we  can  estimate  the  amount  of  model  produced  by  a  unit-valued 
Q'1  at  the  geographic  point  of  interest.  This  technique  is  called  as  point  spreading 
function  (psf)  and  the  rest  of  the  model  other  than  unity  is  assigned  to  zero.  The  synthetic  data  is 
used  as  an  input  model  and  inverted  via  backprojection  to  estimate  the  spreading  of  a  unit  grid. 

For  the  analysis  of  lateral  variation  in  Tibet,  we  first  found  a  profile  which  best  fits  the  station 
points  using  least  squares.  Then,  we  transformed  the  spherical  coordinates  to  Cartesian 
coordinates  and  rotated  using  the  Euler  pole  rotation.  The  purpose  is  to  eliminate  the  errors 
caused  by  the  sphericity  of  the  earth.  The  new  station  coordinate  arranged  to  lie  along  the 
equatorial  latitude.  The  two-station  pairs  for  each  array  (INDEPTH-II  and  INDEPTH-III)  were 
arranged  across  the  profile  so  that  they  do  not  oversample  one  point  and  moving  with  the  same 
separation.  The  profiles  were  inverted  using  a  0. 1  degree  cell  size. 

7.7  RESULTS. 

The  INDEPTH  II  profile  is  located  in  Southern  Tibet  where  partially  molten  middle  crust  exists 
north  of  the  Kangmar  dome,  about  80  km  south  of  the  IYS,  to  central  Tibet  (e.g.  Nelson  et  al., 
1996;  Kind  et  al.,  1996).  We  used  both  Spectral  Inversion  and  Two-Station  methods  allowing 
comparing  the  results  obtained  from  different  estimates  of  Lg  Qo. 
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Figure  7-4.  Example  of  displacement  spectra  with  the  noise  window  from  INDEPTH  II.  Event  1994  236  21:08 
Mb=5.5.  Distance  to  first  station  BB05  is  about  350  km. 

To  perform  the  spectral  inversion  on  the  stations  in  Figure  7-1,  we  first  calculated  the 
displacement  spectra  of  each  station  and  removed  the  effect  of  noise  and  geometrical  spreading 
(Figures  7-2  and  7-4).  Spectra  with  a  good  S/N  ratio  are  selected  for  inversion.  The  frequency 
range  in  the  inversion  is  adjusted  to  minimize  the  noise  contribution.  The  selected  spectra  are 
used  for  the  inversion  after  a  couple  of  successive  iterations  that  are  narrowed  around  the 
solution.  We  select  the  solution  with  the  minimum  residual  (Equation  7.5).  Since  our  main 
concern  is  the  attenuation  of  Lg,  we  did  not  include  M0  and  fcon  the  map,  though  they  are  given 
in  Table  7-1.  Figure  7-5  shows  the  result  of  Lg  Qo  at  Southern  Tibet. 

It  is  notable  to  see  the  higher  attenuation  at  the  north  of  the  Kangmar  dome.  Because  Sn 
propagation  is  efficient  from  nearby  earthquakes  originated  in  the  Indian  Shield  and  the 
Himalayas  (Ni  and  Barazangi,  1983;  McNamara,  1996)  and  the  coda  of  Sn  will  contaminate  the 
Lg  phase,  thus  these  nearby  southern  events  are  not  used  in  this  study.  The  significant  difference 
for  Lg  Qo  is  for  the  event  occurred  in  Indian  Shield  with  Lg  Qo  =225.  Such  a  low  Q  value  is 
clearly  not  a  true  representation  of  the  Indian  Shield;  rather  it  is  a  measure  of  attenuation  through 
the  Indian  Shield,  Himalayas  and  Southern  Tibet. 
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Table  7-1.  Spectral  inversion  results  of  INDEPTH  II. 


Date(yyjdy) 

O.Time 

Qo 

r|(eta) 

Mo(dyne-cm) 

fc(hz) 

94160 

22:53 

136 

0.7 

1.2xl023 

1.3 

94176 

08:32 

110 

0.8 

9  xlO22 

0.6 

94180 

18:22 

112 

0.5 

7  x  1 024 

0.4 

94181 

00:48 

86 

0.7 

1.4  xlO24 

0.4 

94204 

20:57 

74 

0.7 

4.3  xlO24 

0.7 

94205 

23:39 

180 

0.5 

1.1  xlO23 

0.9 

94216 

05:08 

70 

0.8 

9  x  1 022 

0.5 

94236 

21:08 

116 

0.3 

1.1  xlO23 

1.3 

94237 

16:47 

90 

0.6 

6  x  1 022 

0.6 

94243 

23:29 

225 

0.7 

2.6  xlO23 

1.3 

Two-station  Lg  Qo  results  are  shown  in  Figure  7-6.  We  used  these  inter-station  Q  values  to 
measure  lateral  changes  along  the  profile.  Since  the  INDEPTH  II  experiment  was  operated 
about  six  months  it  is  hard  to  find  many  events  aligned  along  the  same  great  circle  path.  We 
used  one  event  (1994  08/24(236)  21:08:47.3  Mb=5.5),  which  occurred  north  of  the  array  (Figure 
7-7).  Although  the  station  separation  should  normally  be  kept  longer  than  1 50-200  km,  we  used 
even  shorter  distances.  The  fact  is  that  dramatically  high  attenuation  allows  us  to  use  short 
station  separation.  The  same  data  preparation  as  for  spectral  inversion  is  applied  to  waveforms 
such  as  removal  of  geometrical  spreading  and  noise  effects.  Each  spectrum  is  then  smoothed 
using  17-point  running-average  (Figure  7-8).  Figure  7-9  shows  the  lateral  variation  of  Lg  Qo 
along  the  INDEPTH  II  profile.  We  observe  very  low  attenuation  areas  just  north  of  the  IYS 
(-50-60).  The  slight  increase  of  Q  (-100)  is  well  observed  for  the  southernmost  part  of  the 
INDEPTH-1I  array  where  is  geographically  known  as  the  high  Tethyan  Himalayas. 


Table  7-2.  Spectral  inversion  results  of  INDEPTH  III. 


Date(yyjdy) 

O.Time 

Qo 

r|(eta) 

Mo(dyne-cni) 

fc(hz) 

98318 

04:11 

ill 

0.5 

2xl022 

1.2 

98331 

06:32 

87 

0.4 

4xl022 

1.7 

98359 

10:44 

86 

0.6 

6xl022 

1.0 

99068 

10:09 

155 

0.3 

1  x  1 023 

0.4 

99092 

09:39 

97 

0.5 

9x  1 022 

0.8 

99149 

20:10 

108 

0.4 

9.4  xlO23 

0.9 
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Figure  7-5.  Spectral  inversion  results  INDEPTH  II.  Red  dots  and  indicate  events,  red  numbers  are  Julian  days.  Lg 
Qo  is  given  in  shaded  area. 
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Figure  7-6.  Station  pairs  calculated  using  17  point  smoothed  spectra.  Since  the  smoothing  procedure  will  effect 
frequency  dependence  term  we  did  not  include  at  results. 
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Figure  7-7.  Location  of  event  1994  236  and  examples  of  waveforms  with  their  spectra.  Note  the  significant 
attenuation/blockage  of  Lg  for  higher  frequencies. 
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Figure  7-8.  The  17  point  smoothed  spectra  of  event  1994  236  21:08. 


INDEPTH  III  recorded  along  a  NNW-SSE  seismic  profile  in  central  Tibet  (Figure  7-1)  across  the 
Banggong-Nujiang  suture  (BNS).  During  the  period  of  one  year  some  regional  earthquakes  were 
recorded.  The  dense  distribution  of  stations  allowed  us  to  use  both  spectral  inversion  and  two- 
station  inversion  methods,  despite  a  number  of  quality  events.  Results  are  shown  in  Table  7-2. 

We  used  all  possible  station  pairs  using  a  station  separation  of  about  200  km  (Figure  7-10).  An 
event  that  occurred  very  close  to  the  BNS  allowed  us  to  visually  compare  the  characteristics  of 
Lg  traveling  to  Southern  and  Northern  Tibet  (Figure  7-11).  No  major  differences  in  the  spectra 
are  observed.  The  two-station  inversion  result  is  shown  in  Figure  7-12.  Apparently,  using  the 
two  events  along  the  same  great  circle  path  for  two-station  inversion  gives  similar  results  with  no 
major  change  along  the  profile.  To  show  the  attenuation  difference  in  central  and  Southern 
Tibet, 
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Figure  7-9.  Lateral  variation  of  Lg  Q0  through  INDEPTH  II  profile.  Green  lines  are  error  bars,  black  dots  are  the 
lateral  resolution  in  units  of  cells. 
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we  selected  two  events  recorded  along  INDEPTH  II  and  INDEPTH  III  arrays.  The  drastic 
change  in  both  amplitude  and  frequency  within  the  same  distance  is  seen  for  INDEPTH  II 
stations  (Figure  7-13).  This  difference  indicates  the  Southernmost  Tibet  where  the  mid-crust  is 
molten  has  a  large  attenuation  on  the  Lg  phase. 

Qo=91.2  ±  13.7  Qo=93.5  ±  23.5  Qo=58.4  ±  1.8  Qo=77.4±8.5  Qo-83.8  ±9.0 


Figure  7-10.  Two  station  pairs  for  INDEPTH  III  event  1998  331. 


For  paths  crossing  the  Tibetan  Plateau,  Lg  phase  is  not  observed  for  paths  longer  than  600-700 
km  (McNamara  et  al.,  1996).  McNamara  et  al.  analyzed  data  of  this  experiment  for  the  presence 
of  Lg  as  well  as  attenuation  to  each  station  using  data  collected  by  the  1991-1992  PASSCAL 
Tibet  experiment.  In  order  to  cover  a  larger  area,  we  included  the  Lg  paths  to  INDEPTH  II 
stations  (Figures  7-14  and  7-15).  Before  selecting  earthquakes  for  Lg  attenuation  we  checked 
the  presence  of  this  phase  at  each  station.  The  stations  and  events  are  shown  in  Figure  7-15, 
events  and  results  are  listed  in  Table  7-3. 
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We  did  not  apply  the  two-station  method  for  detailed  lateral  variation  along  the  array  because  the 
station  distribution  is  so  sparse  for  such  an  application.  However,  Xie  (2002)  has  used  two- 
station  pairs  by  averaging  them  for  a  small  region  and  found  Lg  Qo  =126  for  a  part  of  Eastern 
Tibetan  plateau. 


ST15 
Dist=405km 


ST27 
Dist=460km 


4+ A^Wk- 


Figure  7-11 


.  Comparison  of  two  waveforms  from  south  and  north  of  suture, 
observed. 


There  is  no  significant  change  is 


49 


80 


85 


90 


95 


Figure  7-12.  Station  and  events  with  Lg  Q0  values  for  INDEPTH  III. 


1998  Nov  27  (331) 
06:32  :05.9 


Time(»ec) 


1994  Aug  24  (236) 
21:08  :47.3 


Figure  7-13.  Comparison  of  attenuation  for  the  same  epicentral  distances  of  INDEPTH  II  and  INDEPTH  III. 

Waveforms  are  bandpass  filtered  between  0.5-5  Hz.  Distance  to  the  event  is  given  in  parenthesis  of 
each  station  name. 
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Figure  7-14.  Two-station  lateral  variation  inversion  results  for  INDEPTH  III.  Green  lines  are  error  bars.  Black 
dots  at  lower  part  of  the  plot  shows  lateral  resolution  in  cells  of  0.1  degrees  distance. 


51 


No.  of  cells  (x  0.1  degrees) 


40 

36 

32 

28 

24 

20 

Figure  7-15.  Lg  propagation  in  Tibet  using  INDEPTH  II  and  1991-92  PASSCAL  Experiment. 
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Figure  7-16.  Stations  and  events  used  for  spectral  inversion  of  1991-1992  PASSCAL  Experiment. 
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Table  7-3.  Spectral  inversion  results  of  1991-1992  Tibet  PASSCAL  Experiment. 


DATE(yyjdy) 

O.Time 

Qo 

A (eta) 

Mo(dyne-cm) 

fc(Hz) 

91199 

13:24 

81 

0.7 

3.6x1 023 

0.6 

91201 

19:02 

107 

0.6 

6  xlO22 

0.7 

91245 

11:05 

94 

0.6 

8  xlO23 

0.7 

91251 

23:54 

83 

0.7 

1  xlO23 

2.4 

91263 

11:16 

110 

0.8 

3  xlO24 

0.3 

91323 

1:04 

84 

0.7 

2.8  xlO23 

0.5 

91336 

19:45 

95 

0.6 

21  xlO23 

0.5 

91349 

15:59 

81 

0.7 

6.2x1 023 

0.6 

91357 

1:58 

89 

0.6 

2.5  xlO23 

0.7 

92104 

3:47 

123 

0.5 

9  xlO22 

0.6 

7.8  CONCLUSIONS. 


The  high  attenuation  in  the  crust  of  Tibet  was  estimated  using  three  different  array  dataset.  The 
Lg  Qo  in  the  eastern  part  of  Tibet  is  found  about  ~100-1 10  similar  to  the  two-station  inversion 
study  of  Xie  (2002)  ( Lg  Qo  =126±9).  In  this  study  we  found  that  there  are  regions  of  extremely 
high  attenuating  Lg  in  the  southernmost  of  Tibet.  It  appears  the  high  attenuation  is  intimately 
related  to  the  mid-crustal  partial  melt.  Fan  and  Lay  (2002)  found  higher  Lg  Q  using  earthquakes 
that  are  more  distant  from  station  WMQ.  They  found  Lg  Qo  =316  for  the  southern  part  of  Tibet. 
Our  Lg  Qo  estimates  for  Southernmost  Tibet,  using  data  collected  within  Tibet,  are  very  low  (60- 
100).  The  lowest  attenuation  is  in  the  vicinity  of  IYS  (~60)  and  higher  for  the  area  south  of 
Kangmar  dome  (~100)  (Figure  7-9).  However,  large  error  due  to  limited  data,  which  sampling 
this  complicated  region,  the  location  of  the  change  of  low  to  high  attenuation  in  Southern  Tibet 
cannot  be  resolved  yet.  Moreover,  any  major  structure  change  at  IYS  that  can  cause  large 
scattering  is  not  resolved.  We  observe  that  INDEPTH  III  Lg  Qo  results  are  also  about  80-110. 
There  is  apparently  no  large  attenuation  change  in  central  to  Southern  Tibet.  Data  in  Northern 
Tibet  is  limited  and  how  the  Lg  attenuation  behaves  is  not  well  known.  Very  low  Q  found  in 
Tibet  would  explain  the  limiting  distance  of  600-700  km  for  observing  Lg  in  Eastern  and  central 
part  of  Tibet.  These  values  of  Lg  Qo  may  offer  an  alternative  explanation  for  the  Lg  ‘blockage’ 
along  the  boundaries  of  the  plateau.  The  low  Lg  Qo  values  are  consistent  with  high  temperature 
and  presence  of  melt  and  fluid  in  the  Tibetan  crust.  The  lowest  Qo  in  Southern  Tibet  is  probably 
caused  by  a  molten  mid-crust.  Attenuation  values  given  above  are  very  low  for  continental  areas 
anywhere  and  can  explain  the  high  attenuation  in  Southernmost  Tibet. 
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